Gravitational waves from the evolution of the /-mode instability in neutron stars 
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We study the dynamical evolution of the gravitational-wave driven instability of the /-mode in 
rapidly rotating relativistic stars. With an approach based on linear perturbation theory we describe 
the evolution of the mode amplitude and follow the trajectory of a newborn neutron star through its 
instability window. The influence on the /-mode instability of the magnetic field and the presence 
of an unstable r-mode is also considered. Two different configurations are studied in more detail; a 
standard N = 1 polytrope with a typical mass and radius and a more extreme polytropic N = 2/3 
model which describes a supramassive neutron star. We study several evolutions with different 
initial rotation rates and temperature and determine the gravitational waves radiated during the 
instability. For reasonable values of the mode saturation amplitude, i.e. with a mode energy of 
about 10 -6 Mqc 2 , the gravitational-wave signal can be detected by the Einstein Telescope detector 
from the Virgo cluster. The magnetic field affects the evolution and then the detectability of the 
gravitational radiation when its strength is higher than 10 12 G, while the effects of an unstable r- 
mode become dominant when this mode reaches the maximum saturation value allowed by non-linear 
mode couplings. However, the relative saturation amplitude of the /- and r-modes must be known 
more accurately in order to provide a definitive answer to this issue. From the thermal evolution 
we find also that the heat generated by shear viscosity during the saturation phase completely 
balances the neutrinos' cooling and prevents the star from entering the regime of mutual friction. 
The evolution time of the instability is therefore longer and the star loses significantly larger amounts 
of angular momentum via gravitational waves. 

PACS numbers: 04.30.Db, 04.40.Dg, 95.30.Sf, 97.10.Sj 



I. INTRODUCTION 

After a supernova explosion a nascent neutron star, if 
rapidly rotating, may develop non-axisymmetric insta- 
bilities and radiate a significant amount of gravitational 
waves pQ. These instabilities may originate from differ- 
ent physical processes and grow on dynamical and secu- 
lar timescales. The dynamical instabilities occur in star 
with high rotation rates. In uniformly rotating stars, the 
typical case of the bar mode instability sets is when the 
kinetic to gravitational potential ratio is T/|W^| > 0.27. 
Secular instabilities are driven by dissipative processes 
and thus develop on longer timescales but at smaller ro- 
tation rates. For instance, the viscosity driven bar mode 
instability typically sets in when T/|W| > 0.14 [2], while 
the gravitational- wave driven I = m = 2 /-mode instabil- 
ity requires roughly T/|W| > 0.13 0], and in more com- 
pact stars T/|W| > 0.07 [j. The conditions are much 
better for higher multipole /-modes which can be gravi- 
tationally unstable at even lower rotation rates [3]- 

In rotating neutron stars a non-axisymmetric mode 
may be driven unstable by gravitational radiation via 
the well-known Chandrasekhar-Friedman-Schutz (CFS) 
mechanism [5HZ]. This instability occurs when a mode 
which is counter-rotating with respect to the star is seen 
co-rotating by an inertial observer. The angular momen- 
tum radiated by gravitational waves induces an increas- 
ingly negative angular momentum of the mode which 
thus becomes unstable [SHZ] • The mode's growth ends 
when some dissipative process suppress the instability. 

The mode's amplification expected during the insta- 



bility may generate a significant gravitational wave sig- 
nal which can be potentially observed by the current 
and next generation of Earth based laser interferome- 
ters. When the gravitational waves will be observed, 
the identification of the oscillation modes from the spec- 
trum will help us to unveil the properties of the dense 
matter at super-nuclear densities and clarify the neutron 
star physics. The analysis of the observed Quasi Peri- 
odic Oscillations in Magnetars already provided impor- 
tant insights (HI 12'. and even more relevant results are 
expected when Astereoseismology will be applied to the 
gravitational- wave signal [TJ1 [2] . 

Among the several modes which can be driven unsta- 
ble, the /- and r-modes are the most important due to 
their relatively short growth timescale and an efficient 
emission of gravitational waves. The r-mode instability 
has attracted more attention so far as it is CFS unsta- 
ble at any rotation rate and may grow in about few tens 
of seconds in rapidly rotating stars. There is in fact an 
extensive literature which studied the r-mode instability 
for various stellar models and considered it as a possi- 
ble candidate for limiting the star's rotation below the 
Kepler velocity (see for instance [HI [16] and references 
therein). 

The /-mode instability occurs instead in rapidly ro- 
tating stars and preferable in more compact objects |17j . 
In Newtonian stellar models with stiff equation of state 
(EoS), the neutral point of the I = m = 2 /-mode is close 
to the maximum rotation rate, while in softer equations 
of state this mode is always stable [T8H20] . The multi- 
poles for which the /-mode becomes unstable for smaller 
rotation rates are the I = m = 3 and 4, which have a large 
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instability window and still a reasonable short growth 
time [3] I18H20) . The conditions are more promising in rel- 
ativistic stellar models as the /-mode is driven unstable 
at smaller rotation rates, and also the quadrupole /-mode 
may have a significant instability window [H [XT] [5T] . 

To date the evolution of the /-mode instability has 
been studied only for the I = m = 2 case with non-linear 
dynamical simulations in Newtonian gravity. These 
works used both ellipsoidal configurations (22] [23] and 
compressible stellar models with uniform [24] and differ- 
ential [25! rotation. In these non-linear simulations, the 
effects of viscosity has been neglected and the bar mode 
instability is driven by a gravitational radiation reaction 
term which has been incorporated in the Newtonian dy- 
namical equations. In particular, the strength of this 
Post-Newtonian term has been artificially increased to 
shorten the instability evolution and make feasible its 
analysis within the simulation time. For typical neutron 
stars parameter these non-linear studies find that the 
gravitational-wave signal emitted during the bar mode 
instability may be detected by Advanced LIGO from a 
source located in the Virgo cluster. More importantly, 
while the star spins down the mode's frequency decreases 
toward the more sensitive frequency band (100 Hz) of the 
gravitational- wave detectors. 

In the present work we study the evolution of the 
gravitational-wave driven /-mode instability with a dif- 
ferent approach. It is based on the linear perturba- 
tion framework developed for the r-mode instability 
in [53 |2Z] • From the evolution equations for the mode 
energy, total angular momentum and temperature we de- 
rive a system of ordinary differential equations to de- 
scribe the mode's amplitude, the star's rotation and the 
thermal evolution. In this analysis we include the effects 
of viscosity, magnetic field, and consider also the impact 
of an unstable r-mode on the /-mode instability. The 
coefficients of these evolution equations depend on the 
stellar model and mode properties. More specifically, for 
a given neutron star's model we have to determine the 
/-mode frequency and eigenf unctions, and with these in- 
formation calculate the relevant viscous and gravitational 
radiation timescales. 

We consider rapidly rotating and relativistic models 
with uniform rotation and extract the /-mode properties 
(frequency and cigcnfunctions) from the time evolutions 
of the linearised dynamical equations. In particular, we 
simplify the problem by using the Cowling approxima- 
tion, i.e. we neglect the space-time perturbations in our 
linearised problem. The accuracy of this approximation 
is to better than 20% for the quadrupole /-mode, but 
increases considerably for higher multipoles. In a sec- 
ond step, the mode's frequency and eigenfunctions are 
inserted in appropriate volume integrals which determine 
the viscous damping times and the gravitational radia- 
tion growth time [IS] . 

We study in more detail two polytropic neutron star 
models with different compactness and maximum rota- 
tion limit. The first is a standard star with polytropic 



index N = 1 and baryonic mass Mb = 1.4M Q . The sec- 
ond is an N = 2/3 polytrope with Mb = 1.6M Q , which 
represents a supramassive star, i.e. a configuration that 
is stabilised by rotation and that in the non-rotating limit 
does not have a static stable configuration [2"5ll30| . 

Considering several configurations we find that the 
gravitational-wave signal generated by an unstable /- 
mode may be potentially detected with the Einstein Tele- 
scope (ET) from a source located in the Virgo cluster. 
For instance, the gravitational characteristic strain gen- 
erated by the I = m = 3 and 4 /-mode is shown in Fig- 
ure [I] for the N = 1 and N = 2/3 polytropic models with 
a relatively weak magnetic field, B p = 10 11 G. This char- 
acteristic strain is determined by integrating in time the 
signal generated by an /-mode during the instability and 
setting a maximum saturation energy of E ~ 10 -6 Mqc 2 . 
However in order to not overestimate the signal when the 
instability evolves on long timescales, we have considered 
a maximum detector integration time of 1 year. For more 
details on these results and the effects of higher magnetic 
fields and r-modes see Sec. [V] 

A general expectation is that superfluidity should fur- 
ther restrict the parameter space of the instability. In 
fact, if a star cools down below the transition tempera- 
ture where the neutrons of the core become superfluid, 
the /-mode instability should be efficiently suppressed by 
the mutual friction force [3T]. However, for typical stel- 
lar parameters we find that the heat generated by shear 
viscosity during the saturation phase of an unstable /- 
mode may counterbalance the neutrino cooling. In fact, 
our results show that an unstable star may follows a quasi 
isothermal trajectory within the instability window with- 
out cooling below the critical temperature of the core's 
neutrons. As a result the evolution time is longer and 
the star loses more angular momentum via the emission 
of gravitational waves. 

We present the formalism for studying the /-mode in- 
stability in Sections [IT] and [In] while in the Appendix we 
provide the equations for the gravitational-wave instabil- 
ity driven simultaneously by two modes. The results for 
the instability evolution are described in Section [IV] and 
the gravitational- wave signal is determined in Section [V] 



The concluding remarks can be found in Section VI 



II. THE /-MODE INSTABILITY 

The standard approach for studying secular instabili- 
ties is based on linear perturbation theory and consists 
of two distinct steps. First, we determine the main prop- 
erties of a mode (frequencies and eigenfunctions) , which 
will be inserted in a second step into appropriate vol- 
ume integrals for the estimation of the relevant timescales 
due to viscosity and gravitational radiation. Since these 
timescales are generally much longer than the oscillation 
period of the mode, their properties can be determined 
by solving the inviscid problem |19] . 

In this work, we consider relativistic models of rapidly 
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FIG. 1: Characteristic strain generated by the /-mode instability for a polytrope with B v = 10 11 G . The source is located 
at 20 Mpc and the saturation amplitude of the /-mode is set to a sat = 10 -4 (see Sections III A and IV A I. Left panel: The 
gravitational- wave signal emitted by the I = m = 4 /-mode for the two stellar models with TV = 1 and TV — 2/3. The number 
near the vertical lines denotes the initial rotation rate Right panel: The signal emitted the by I = m = 3 /-mode for the 

more massive model with TV — 2/3. The sensitivity curves of Advanced LIGO and the Einstein Telescope (ET) are shown in 
both panels. The gravitational-wave signal generated during the /-mode instability may be detected by ET for the most part 
of the instability window of the TV = 2/3 model. 



rotating neutron stars and determine the /-mode fre- 
quency and eigenfunctions by post-processing time evo- 
lution data of the relativistic perturbation equations. All 
the details of the formalism and numerical methods for 
studying the oscillations of rapidly rotating stars is given 
in [31I33]. 

In a realistic star, an oscillation mode approximately 
evolves as e lwt ~*/ T , where oj is the mode's frequency and 
r the damping/growth timescale. A way to determine r 
is to calculate the variation in time of the mode's energy 
given in the rotating frame by 



E = 



1 



dV 



pSu l 5u* + — (Spdh* + Sp*Sh) 



(1) 



where the asterisk denotes a complex conjugate and the 
integral extends over the volume of the star. The quan- 
tity p represents the background mass density, while 5p, 
6h and 8u l are perturbations of the mass density, the en- 
thalpy and the velocity, respectively. The latin indices 
denote the spatial components of a generic four vector. 
Equation ([!]) is also defined as the canonical energy in 
the rotating frame [54"] . 

As long as the energy of a mode depends quadratically 
on its perturbations, one can show that 



dE 

~dt 



2E 

T 



(2) 



where the global timescale r of equation ^ results from 
the individual dissipative mechanisms that act on the 



mode. For stars in which only gravitational waves, shear 
and bulk viscosity dissipate energy, we have [19j : 



1 



' gw 



1 



1 

Tb 



(3) 



The gravitational radiation timescale can be calculated 
by using the standard multipole formula 1 [35] . 



1 



LO 

2E 



Y,Ndoj- mn) 2l+1 (\SD ln \ 2 + \SJ b 



l>2 



(4) 



where uj is the mode frequency as measured in the rotat- 
ing frame, 5Di m , SJi m are the mass and current multipole 
moments and TV/ is given by 



TV, 



47rG (I + + 2) 
c^+^Z (Z - 1) [(2Z + l)!!] s 



(5) 



Since pressure modes mainly radiate via the mass multi- 
pole moments which are determined by 



SD, 



dVr l 8pY; ml 



(6) 



1 In our calculations we correct the r gw determined via equa- 
tion |4jl with a constant factor. This procedure provides a bet- 
ter agreement with the relativistic results (no Cowling approx- 
imation) of the gravitational radiation timescale; see Doneva & 
Kokkotas (in preparation) for more details. 
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the current multipole moments are neglected for comput- 
ing damping times of the /-mode. 

The dissipative timescales due to bulk- and shear- 
viscosity can be determined by the following volume in- 
tegrals 

~ = ~/^C^*, (7) 
^^JdVvS^Sa^, (8) 

where C and r\ are the bulk and shear viscosity coeffi- 
cients, while the stress tensor Sa^ is given in terms of 
velocity perturbations, 

Sa^ = * (Vw + V'Su* - ~g id VSa\ , (9) 
8a = Vj5u j , (10) 

and gij is the spatial metric tensor. 

The bulk and shear viscosity coefficients depend on 
the state and composition of the neutron star matter. 
For a neutron, proton, electron (npe) matter in normal 
state, i.e. with no superfluid/superconducting compo- 
nents, these coefficients can be written in a simple ana- 
lytical form (T9l [36] . More specifically, the bulk viscosity 
coefficient is given by [37] 

C = 6 x 1(T 5 V T 6 g cm-V 1 , (11) 

where T is the star's temperature. This expression is 
strictly valid for "small" oscillation amplitudes, within 
the so-called sub-thermal regime. If the mode ampli- 
tude is sufficiently large, non-linear terms may become 
important and will increase the strength of the bulk vis- 
cosity 38, 39 . However, the impact of the non- linear 
bulk viscosity on the /-mode instability is moderate and 
appears at rather large mode amplitudes [40]. In the 
present work we limit the /-mode growth to compara- 
tively small amplitudes and then neglect the effects of 
non-linear bulk viscosity. 

In normal npe matter, shear viscosity is dominated 
by neutron collisions and can be parametrized with the 
following coefficient [It?! [55] : 

?7 = 347p 9/4 T- 2 gcm- 1 s- 1 . (12) 

For an inviscid star, there is now the following situa- 
tion: if the star is oscillating at a sufficiently high rota- 
tion rate, the condition for the secular CFS-instability, 
bj (lj — mil) < 0, might be fulfilled. In this case, the 
stellar rotation drags a counter-rotating mode into coro- 
tation as seen from an inertial observer, leading to a grav- 
itational wave driven, exponential grow of the initial per- 
turbation (it is r gw < then; see equation Q). However, 
the other dissipative mechanisms described above also 
operate in realistic stars and tend to stabilise the mode 
evolution. The instability onset and the mode growth is 
therefore controlled by the global timescale which has to 



be t < 0. This global timescale is generally a function of 
the star's rotation rate f2 and temperature T so that the 
instability has a natural representation in a T — plane, 
where the region above the critical curve for r = de- 
notes the so-called instability window (e.g. see Figure|2|. 

This description of the instability is clearly static as 
no information about the evolution of the star and mode 
amplitude is provided. The evolution problem will be 
addressed in the next Section. 



III. EVOLUTION EQUATIONS FOR THE 
/-MODE INSTABILITY 

To study the evolution of the /-mode instability we de- 
rive a system of equations by using and extending the for- 
malism already developed for the r-mode instability [26) . 
Our equations consider also the effects of thermal heat- 
ing, magnetic torque and the impact of an unstable r- 
mode on the evolution of the /-mode instability and its 
gravitational-wave signal. 

In the non-linear study of the secular bar-mode insta- 
bility [22TI25] . it is reasonable to discern two dynamical 
phases of the /-mode instability, namely the mode's ex- 
ponential growth and its non- linear saturation [26] . As 
the star enters the instability window, the mode grows 
exponentially while the star slowly spins down on viscous 
timescales. This initial growth phase ends either when 
non-linear dynamics or dissipative processes saturate the 
mode amplitude at a finite value. After this point, the 
mode amplitude is nearly constant while the star loses 
angular momentum via gravitational radiation [2H [35] . 

The dynamical equations that describe these two 
regimes can be derived using the relations for the canon- 
ical mode energy (in the rotating frame) and the stellar 
angular momentum. In this section we focus only on the 
/-mode instability, while in the Appendix we also address 
the more general case in which two modes can be simul- 
taneously driven unstable by gravitational radiation. 



A. Mode evolution and star's spin-down 

The mode energy can be considered as a function of 
the mode amplitude and the stellar rotation rate. Hence, 
we may assume that E = aE (f2), where the parameter 
a denotes the mode amplitude and E (f2) describes the 
angular velocity dependence. It is worth noticing that in 
contrast to [55J [30] we define a from the mode energy, 
hence the fluid perturbations scales as ~ a 1 / 2 . After 
taking the time derivative, equation ([2| can be written 
as 

da d In E dQ 2a 

M +a ^T^ = -V (13) 

A second equation can be derived from the evolution 
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of the angular momentum, 



and the coefficients are given by the following expres- 



dJ 
~dt 



dJ a 



dJ„ 



dt 



dt 



(14) 



where J = J s + J c is the total angular momentum, 
which consists of the star's angular momentum J s and 
the canonical angular momentum of the mode J c . The 
latter one can be related to the mode energy by the well 
known equation |34j : 



Jc 



-E. 



(15) 



The dissipative terms in equation ( 14 1 are the gravita- 
tional radiation torque, 



dJ a 



dt 



2J ( 

T, 



(16) 



gw 



and the magnetic torque dJ mag /dt. 

We consider a standard oblique rotator model in vac- 
uum with a dipolar magnetic field, which has [4Tj 



dJ„ 



■ 2 

sin x 



dt 



6c3 



(17) 



where R is the stellar radius, B p is the magnetic field at 
the magnetic pole and \ is the inclination angle between 
the rotation and magnetic axes. For simplicity we assume 
an orthogonal rotator (\ — 7T /2). This configuration pro- 
vides also a magnetic spin-down comparable to that of a 
standard pulsar model with magnetosphere [HI [42] . 



From equation (15) and the functional form of E we 



can write the canonical angular momentum as J c = 
aJ c (O). In this way we can simplify the calculation and 
determine from equation (14) the following expression: 



da 1 / dJ s dJ c \ dfl ( a 1 

dt J r \d£l dfl J dt \ r gw T ma 

where we have introduced a timescale for the magnetic 
torque, 



' mag 



1 dJ ma g 

2J C dt 



(19) 



From equations ( 13 ) and ( 18 ) we can now derive a 

(20) 



system of ordinary differential equations: 

da 2a 2a 1 + aQ 2P a 



dt T gw t v D 

dQ. _2F ( a 1 



D r„ 



where the total viscous damping time is defined as 



1 _ 1 1 

T v T s T b 



(21) 



(22) 



_ dJc ( dJs_ 
Q ~~ ~dtt \dtt 

d\nE 
P = F. 

dn 



F = 3 >{W ■ < 23 > 

D = 1 + a (Q - P) . (24) 



The form of equations (20) and (21 ) agrees with the evo- 



lution expected during the initial phase of the instability, 
in which the /-mode grows exponentially as a ~ e*^' 1 " 9 ™' 
and the star slows down on viscous timescales. 

As we have already pointed out, we expect that dur- 
ing the non-linear saturation regime the amplitude of the 
mode remains nearly constant [23J [5S]. We can there- 
fore approximate this evolution phase by assuming that 
da/dt = 0, which allows us to recast equations ( 20 1-( 21 ) 
into a single relation, 



dn 

lit 



2F 



1 + aQ \ t. 



1 



• gw 



' mag 



(25) 



This expression shows that for neutron stars with a small 
magnetic field the spin down during the non-linear satu- 
ration phase is dominated, as expected, by gravitational 
radiation. 



Equations (20)-(21) and (25) describe the time varia- 



tion of the mode amplitude and the star's angular veloc- 
ity, but do not provide any information about the ther- 
mal evolution. In the next section, we address the cooling 
problem and derive the equations for the thermal balance 
of the star. 



B. Thermal evolution 

Few minutes after a core collapse, a neutron star be- 
comes transparent to neutrinos which dominate the cool- 
ing for at least 1000 years [41] . However, an isolated 
neutron star may be re-heated by viscosity and magnetic 
field decay in later stages of its life. The evolution of the 
magnetic field is not well known in stars younger than 
10 4 yr [3S], and it is typically described by approximate 
relations [4"5"H4"T] . For instance, an equation used in [4"5] 
reads 



B = B a 1 



t 

Tde, 



(26) 



ec the field 
10 4 



where B is the initial magnetic field and r^. 
decay timescale. This is approximately Td ec ~ ttn yr 
and it is much longer than the typical evolution of the 
/-mode instability. Hence, the magnetic field decay can 
be reasonably neglected during the unstable phase of the 
/-mode. 

The two viscous processes that operate in our models 
have an opposite effect on the thermal evolution. Bulk 
viscosity dissipates the mode energy by neutrino emis- 
sion, as a result of the Urea reactions that attempt to 
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FIG. 2: Evolution of the I = m = 4 /-mode instability for a relativistic polytrope with N = 1, a baryonic mass of Mb = 1.4 Mq 
and i? p = 10 11 G. The star enters the instability window at different spins (see right panel) and correspondingly the left panel 
depicts the time-evolution of the stellar rotation rate (top left panel) and mode amplitude (bottom left panel), where the initial 
amplitude a — 10 -6 saturates at a sat — 10~ 4 . The shaded region in the right panel represents the temperature range where 
the neutrons and protons of the core are expected to be in a superfluid/superconducting state in accordance with recent models 
for the observed cooling of Cassiopeia A [431 144] . 



restore the beta equilibrium in a displaced fluid element. 
These neutrinos induced by mode oscillations escape at 
infinity and cool down the star. However, this effect is 
very small compared to the dominant cooling rate of the 
background star and can be therefore neglected. More 
important instead is shear viscosity which dissipates en- 
ergy through particle collisions. The heat generated by 
shear viscosity is proportional to the mode amplitude and 
it can be relevant during a mode instability. In fact as 
we will show in Section |IV| shear viscosity may reduce 
and even completely balance the neutrino cooling of an 
unstable star. 

The thermal evolution of a neutron star can be stud- 
ied with a global energy balance between the relevant 
radiative and viscous processes [4"T] : 



C dT - L 
at 



(27) 



where C v is the total heat capacity at constant volume, 
L v is the neutrino's luminosity produced by the modi- 
fied Urea processes that operate in the background star, 
and H s is the heating rate generated by shear viscosity. 
Although equation p7| ) is strictly valid for isothermal 
neutron stars we consider it a reasonable approximation 
for this work. 

The heat capacity depends in general on the neutron 
star EoS and temperature. An expression frequently used 
in literature is given by |41j 



C v = 1.2 x 10 



:','.) 



M 



P 



-2/3 



10 9 K 



erg K- 1 , (28) 



which has been strictly derived for an ideal Fermi gas 
of neutrons in which p < Ip^, where po is the nuclear 
saturation density. 

The neutrinos' luminosity can be determined by the 
following equation [IT] : 



L v = 5.3 x 10 



:if> 



M 



P 



-1/3 



(— V 

I io 9 k; 



erg s 



m q \ Po j V!o 9 k; ' (29) 

which has already been used for the r-mode instabil 
ity [26] , 



The last term we need to specify in equation (271 is 



the heating rate induced by shear viscosity, which can be 
described by the following equation [2"rj] : 



H„ = 



2E 



(30) 



We now have all the ingredients to study the evolu- 
tion of the /-mode instability with the system of equa- 
tions pOl-fpll) and p7|. 



IV. RESULTS 

In this Section, we discuss the dynamical evolution of 
the /-mode with and without the impact of an additional 
magnetic field. Furthermore, we study also the effect 
of an unstable r-mode on the evolution of the /-mode 
instability. For this, we model the compact objects as 
relativistic neutron stars which obey a polytropic EoS, 

p = Kp 1+1 ' N (31) 
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where p is the fluid pressure, K is the polytropic constant 
and N is the polytropic index. The rest-mass density p 
is related to the fluid energy density e via the relation 
e = p + Np. 

To study the dependence of the /-mode instability on 
the neutron star model, we consider two sequences of 
uniformly rotating stars up to the mass shedding limit 
fix- The first model represents a "standard" neutron 
star with mass Mf, = 1.4 M and polytropic index N = 1. 
The non-rotating member of this sequence has a circum- 
ferential radius of R = 14.15 km, while the maximum 
rotation rate is vk — ^k/27t = 673.1 Hz. 

The second model is more extreme and describes a 
supramassive star, i.e. an object that does not have a sta- 
ble configuration in the static limit [2"5H3T)] . These stars 
with high masses are stabilised by rotation and might 
originate from a gravitational collapse of a massive stel- 
lar core. When a supramassive star spins down it will 
reach a point beyond which gravitational attraction pre- 
vails over the centrifugal force and triggers a gravitational 
collapse to a black hole. We model a supramassive star 
with a polytrope with N = 2/3 and mass Mf, = 1.6 Mq. 
The respective Kepler limit is vk = Ok/27t = 1783 Hz. 

In order to extract physical results from our simula- 
tions we must normalise the mode amplitude a. In this 
work, we use the following definition 



E = aE IC 



(32) 



where -Erot is the stellar rotational energy. For our poly- 
tropic models ana = l corresponds to E ~ 10~ 2 Mqc 2 . 
The definition ( 32 1 implies that the fluid variables and 
the gravitational- wave strain scale as a 1 / 2 . 

The numerical scheme used to study the evolution of 
the /-mode instability is based on a fourth order Runge 



Kutta algorithm. The numerical code evolves separately 
the two phases of the instability. After having specified 
the initial conditions for a, £1 and T, we evolve the initial 



growth phase of the /-mode according to equations ( 20 1- 
( 21 ). When the mode amplitude reaches a predetermined 
saturation value a sa t, the non-linear saturation phase is 

until the star leaves the in- 

In 



(251 



ic numerical simulation ends. 



studied with relation 
stability window and t 
both phases of the instability we use equation (27 1 for 
the thermal evolution. 



A. The evolution of the /-mode instability 

We first consider the evolution of the /-mode instabil- 
ity in a low magnetised star (B p < 10 11 G), and study the 
effect of the magnetic torque and r-mode in the following 
two Sections. 

The ideal conditions for the /-mode instability to work 
are high temperatures and rotation rates. At birth, the 
temperature of a neutron star is around T ~ 10 11 K and 
drops down rapidly in the following few days. We there- 
fore expect that as a rapidly rotating star cools down, 
it enters the /-mode instability window from the high 
temperature side (see Figures [2] and [3). In our numerical 
scheme, we prescribe this behaviour by choosing appro- 
priate initial conditions for VL and T. In particular, for a 
given rotation rate Q we specify the highest temperature 
allowed by the /-mode instability (e.g. see again Fig- 
ures [2] and [3| . The initial mode amplitude is set to a — 
10~ , this corresponds to an energy of E ~ 10 -8 M c 2 
which is a typical value for the total energy loss in grav- 
itational waves due to quadrupole deformations shortly 
after a gravitational core collapse [451 1^5] . Although the 




initial pulsation energy of the I = m — 3 and 4 /-modes 
might be smaller, our results do not change significantly 
for different initial values of a, as during the initial phase 
of the instability the mode's amplitude grows exponen- 
tially. 

Another parameter one needs to specify in the evolu- 
tion scheme is the saturation amplitude of the /-mode. 
From non-linear dynamical simulations one can deter- 
mine upper limits on the quadrupole /-mode [50] , but the 
problem is still open for higher multipole /-modes. We 
then consider the saturation amplitude as a free param- 
eter. In our simulations we choose a sa t — 10~ 4 , which is 
equivalent to E ~ 10~ 6 M Q c 2 . These values are "reason- 
able" and much smaller than what is frequently used in 
literature [23] [26] , where a sa t = 1 • 

For the N = 1 model we consider four cases for the 
I = m = A /-mode instability with different initial ro- 
tation rates and temperatures. Figure [2] shows the re- 
sults for the star's spin, the mode's amplitude and the 
star's trajectory through the instability window. The 
entire evolution lasts for about 500 years, while the du- 
ration of the initial growth phase depends on the ini- 
tial rotation rate. For a star that becomes unstable 
at f2 = SIk (fi = 0.97 Ok) the growth phase lasts for 
about 1 day (6 years). This phase is dominated by the 
growth time r gw which increases gradually for slower ro- 
tating stars, from r gw ~ 10 4 s at the mass shedding 
limit to r gw ~ 10 8 — 10 9 s at the bottom of the insta- 
bility window, SI ~ 0.95 SIk [ITj- After the mode sat- 
urates, Figure [2] (right panel) shows that the star spins 
down at nearly constant temperature (T ~ 10 9 K) as 
a result of the heat generated by shear viscosity which 
completely balances the neutrino cooling. This effect 
prevents the star from entering the superfiuid transition 



zone and increases the duration of the instability. In fact, 
for T < T cn the neutrons of the core become superfiuid 
and mutual friction damps the /-mode very efficiently on 
short timescales [3T]. For the superfiuid critical temper- 
ature we chose a value of T cn ~ 5 x I0 8 K which has 
recently been determined from the cooling curves of Cas- 
siopeia [431144]. 

In Figure [2] (right panel) it is noticeable that the tra- 
jectory of a star which becomes unstable at SI = 0.97 SIk 
is slightly different from the other cases. This is due to 
the slower growth phase of the mode which delays the 
effect of the shear viscosity heating on the thermal evo- 
lution. Therefore, the star first cools down to a minimum 
temperature, then is re-heated by shear viscosity and fi- 
nally evolves quasi-isothermally in the last part of the 
trajectory. 

The evolution of the I = m = 4 /-mode instability for 
the N = 2/3 model is qualitatively similar to the N = 1 
case, but it develops faster and within a significantly 
larger instability window (see right panel of Fig. [3| . At 
the mass shedding limit the gravitational growth time 
is only r gw ~ 250 s in contrast to the r gw ~ I0 4 s for 
the N = 1 model, and the total evolution now lasts 
for about 40 years. The exponential growth of the /- 
mode takes about 5 minutes (14 days) for a star with an 
initial SI = SIk (SI = 0.93 SIk)- Furthermore, Figure [3] 
shows that when a star becomes unstable at lower rota- 
tion rates near the minimum of the instability window, 
e.g. SI ~ 0.89 S^k, the mode amplitude remains small and 
the shear viscosity heating does not balance the neutron 
star's cooling. 

For the N = 2/3 model we have a relevant instability 
window also for the I = m = 3 /-mode. In fact, at the 
Kepler limit the gravitational growth time, r gw ~ 135 s, 
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is even shorter than the I — m — 4 /-mode but increases 
more rapidly with decreasing rotation. The resulting in- 
stability region is shown in Figure [4] for a star that be- 
comes unstable at different rotation rates. The total in- 
stability evolution time is about 30 years and the dura- 
tion of the growth phase is comparable to the I = m = 4 
/-mode case. 

The properties of the instability window and gravita- 
tional growth time make the N = 2/3 massive star a 
better gravitational-wave source than the N = 1 model 
(see Section |v| . 



B. Magnetic torque 

The main mechanism that gradually slows down a neu- 
tron star is electromagnetic radiation, which is powered 
by the star's rotational energy. Since the magnetic torque 
increases considerably with rotation, as shown by the 
dipole formula (17), the magnetic braking may affect the 



evolution of the /-mode instability by accelerating the 
spin-down. This means that an unstable /-mode may 
not have the time to grow significantly and generate a 
detectable gravitational-wave signal. 

We study the impact of the magnetic field on the I — 
m == 4 /-mode instability for the N = 1 and N = 2/3 
models. For each model we consider two dipolar magnetic 
field configurations with, respectively, B p = 10 12 G and 
10 13 G, and focus on an orthogonal rotator (see 



Section III A ) 



We evolve various initial conditions with different tem- 
peratures and rotation rates. The evolution of the mode 
amplitude and the star spin-down is shown in Figure [5]for 
the N = 2/3 model. For a magnetic field of B p = 10 12 
G the total duration of the instability is shorter than 
a factor of about 50 with respect to the B p = 10 11 G 
model. In fact, the evolution lasts about 1 yr (10 yr) 
for the N = 2/3 [N — 1) model. In slower rotat- 
ing models, the magnetic torque clearly dominates the 
evolution even during the initial phase of the instability 
and limits considerably the growth of the mode ampli- 
tude. This behaviour is for instance evident in Figure [5] 
for a N — 2/3 model with an initial rotation rate of 

= 0.89 J1k- The magnetic torque dominates completely 
the evolution when B p — 10 13 G. The unstable star is 
quickly spun down and the mode amplitude is strongly 
limited even in the more massive model (see Figure pi) . 



C. The /-mode versus the r-mode instability 

In rapidly rotating neutron stars several modes can be 
driven unstable by gravitational radiation at the same 
time. The most important ones are definitely the /- 
and r-modes, as they have a comparatively short growth 
time and therefore can generate a significant gravita- 
tional wave signal. A relevant difference between these 
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FIG. 5: The impact of the magnetic torque on the evolution 
of the I = m = 4 /-mode instability for the N = 2/3 model. 
The magnetic field is dipolar and orthogonal to the rotation 
axis. The two top panels display the mode amplitude (left 
panel) and the star's evolution through the instability region 
(right panel) for B p — 10 12 G. The same quantities are de- 
picted in the two lower panels for B p — 10 13 G. 



two classes of modes is that the /-mode only gets un- 
stable in very rapidly rotating stars while the r-mode 
is CFS unstable at any rotation rate. In addition, the 
growth time is typically shorter for the r-mode which 
consequently has a larger instability window. It is then 
reasonable to think that the r-mode should dominate 
the evolution of the gravitational wave driven instabil- 
ity. However before drawing any secure conclusion, it 
is necessary to know the maximum amplitude that each 
mode can reach during the instability. 

Non-linear perturbation calculations, based on the am- 
plitude formalism, show that the r-mode may transfer 
energy to other inertial modes through non-linear mode 
coupling. As a result the maximum amplitude of the 
r-mode may be limited to c a = 10~ 3 — 10~ 5 , where 
E = c? a E mt |5TJ E2- In our notation, these values 
correspond to a saturation amplitude of a r — c 2 = 
10~ 6 — 10~ 10 , since we parametrize the mode energy ac- 
cording to E — a r E Iot . 

As it was already pointed out, the maximum ampli- 
tude of the /-mode is still uncertain. This forces us to 
consider a "reasonable" value for the /-mode saturation 
which we set to a sa t = 10~ 4 ; the same value that was 
used in the previous sections. We then study the com- 
bined evolution of /- and r-mode for different r-mode 
saturation amplitudes. In order to address this problem 
we derive an additional set of equations which are given 
in the Appendix. 

We focus on the I = rn = 2 r-mode and the I — m — 4 
/-mode, which are the most relevant unstable modes in 
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FIG. 6: The effects of an unstable I — m = 2 r-mode on the evolution of the I = m = 4 /-mode instability for a relativistic 
polytrope with N = 2/3, Mj = 1.6 Mq and B p — 10 11 G. The saturation amplitude of the r-mode is varied between a r = 
10" 9 - 10" 5 while the maximum amplitude of the /-mode is kept constant at a sa t = 10 4 . The left panel shows the evolution 
of the stellar rotation (top panel) and the /-mode amplitude (lower panel) for the various examples. The filled-circle in the 
lower left panel denotes the end of the f-mode instability for a r = 1CP 7 . The right panel shows the evolution through the 
instability window, where the dashed line denotes the high temperature part of the r-mode critical curve (r r = 0). 



their respective classes [2] and study the more massive 
stellar model with N — 2/3 and a low magnetic field 
of B p = 10 11 G, as it is potentially a better source 
for gravitational-wave detection. The properties of the 



/-mode have been already determined in Section IV A 



while for the r-mode one needs to consider some approx- 
imations of its viscous damping and gravitational radi- 
ation growth times. More precisely, the various dissi- 
pative timescales are determined via the analytical rela- 
tions given in |53] which have been derived for a uniform 
density star. For the model under consideration here, 
the I = m = 2 r-mode growth time at the mass shed- 
ding limit is r gw ~ 1.3 s, which is two orders of magni- 
tude smaller than the I = m = 4 /-mode growth time 
(fg W = 244.7 s). The oscillation frequency in the rotat- 
ing frame of the I = m = 2 r-mode is determined at 
leading order by uj — 2f2/3, and for simplicity it is also 
assumed that the canonical angular momenta of r- and 
/-mode are equal. These approximations prevent us from 
determining the properties of the r-mode from numerical 
evolutions. However, since the r-mode grows much faster 
than the /-mode it is not expected that the qualitative 
behaviour of the results will change considerably with a 
more accurate description of the r-mode. 

As the star cools down, it is expected that the r-mode 
gets unstable before the /-mode as its instability win- 
dow extends towards higher temperatures (see Figure [6]) . 
The r-mode then quickly reaches its non-liner saturation 
value a s r at and spins down the star. However if the r- 
mode growth is limited by non-linear mode coupling, the 
star evolves towards the critical curve of the /-mode and 



eventually also this mode is driven unstable by gravita- 
tional radiation. 

In Figure [6] we show a representative case in which 
the star enters the critical curve of the r-mode at Q, — 
0.98 f2j< with an initial amplitude a r = 10~ 10 . We con- 
sider several simulations where the /-mode grows up to 
a maximum amplitude of a s r at = 10~ 9 — 10 -5 , while for 
the /-mode we set a sa t = 10~ 4 . The impact of the r- 
mode on the star's spin-down and the total evolution of 
the /-mode is evident already at low saturation a r . The 
star leaves the instability window of the /-mode after 
an evolution time given by t ev ~ 20 x (lO -9 /cv™*) yr. 
When a r > 10 -5 the evolution is completely dominated 
by the r-mode and the /-mode amplitude is strongly con- 
strained to small values. However, for a r < 10~ 7 the 
/-mode still has time to grow and potentially generates 
a detectable gravitational- wave signal (see Section |v| ). 
As expected, it is therefore crucial to know more accu- 
rately the relative saturation amplitude between these 
two modes. This is an interesting aspect that must be 
clarified in a future work. 



V. GRAVITATIONAL WAVES 

The instability of the /-mode is potentially a strong 
source of gravitational radiation and it is then logical 
to estimate the detectability prospectives by the current 
and next generation of Earth-based laser interferometers. 

The two independent polarization states of the grav- 
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FIG. 7: The characteristic strain generated by the /-mode instability for a polytropic star with B p = 10 12 G. The left panel 
shows the gravitational- wave signal of the I — m — 4 /-mode for the N — 1 and TV = 2/3 models, while the right panel depicts 
the signal of the I — m — 3 /-mode for the TV = 2/3 model. The source is located at 20 Mpc and the saturation amplitude of 
the /-mode is set to a aat = 10 -4 . The notation used in this Figure is the same as in Figure [l] 



itational wave strain can be expressed in terms of the 
spin-weighted spherical harmonic _ 2 Y lm as 



h l SS - ih l Z = h lm ,Y l 



where h is defined as 



nl + 2 



(33) 



(34) 



and u>i = uj — mQ is the mode frequency in the inertial 
frame. The coefficient /?; is given by 



8tt 



'(1 + 2)1 



I (I - 1) (21 + 1)!! V (i-2)! 



(35) 



A more suitable quantity to evaluate the gravitational- 
wave detection is the characteristic strain, which also 
takes into account the statistical amplification of the sig- 
nal due to the number of oscillations accumulated in a 
given frequency bandwidth. The characteristic strain is 
defined by the following equation: 



h lm = (\h lm _ 3 Y lm \)y/N t 



eye 



(36) 



where (...) denotes an averaging over the angles (0,(f>). 
The number of oscillation cycles can be expressed as 
N C yc = v T~ev , which is written in terms of the mode 
frequency v = ujj/2tt and the time spent near a given 
frequency r ev . 

For magnetic fields with a magnitude of B p < 10 11 
G, the impact of the magnetic torque on the /-mode 
instability is irrelevant (see Section IV I, and the mode 



evolves like in an unmagnetised star. This means that 



during the initial exponential growth of the mode the 
star spins down on viscous timescales, i.e. is nearly 
constant, and the radiated gravitational- wave signal is 
virtually monochromatic. As a result, the quantity r ev 
can be approximated by the accumulated evolution time 
of the star during this initial phase with Q ~ const. Af- 
ter the /-mode saturates, the amplitude remains approx- 
imately constant and gravitational radiation removes an- 
gular momentum from the star. This leads to a frequency 
variation which can be well approximated by the follow- 
ing relation: 



dt 
dp 



(37) 



When a star has higher magnetic field strengths, B p > 
10 11 G, the magnetic torque affects also the initial phase 
of the instability and accelerates the spin-down. The /- 
mode frequency therefore evolves as the star loses angular 
momentum by magnetic braking and equation (37) has 



to be used throughout the entire instability evolution. 

However, realistic calculations of the characteristic 
strain also have to consider the maximum signal inte- 
gration time allowed by the detector technology. We 
assume that a detector may integrate the signal at most 
for 1 year, and consequently calculate h l ™ for r ev < lyr. 
This means that in N cyc we use the evolution time of 
our simulations whenever r ev < lyr, otherwise we set 
T ev = lyr. 

The results for the N — 1 and N = 2/3 models are 
shown in Figure [T] and in Figures [7J [8] together with the 
sensitivity curves of Advanced LIGO and ET. The differ- 
ent noise curves of the detectors are determined by using 
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hrms = \J vSh{v), where Sh{v) is the power spectrum of 
the detector at hand. 

For each model we determine the characteristic strain 
of the various evolutionary paths studied in Sec. |IV| in 
which the onset of the /-mode instability happens at dif- 
ferent rotation rates. As shown in Figure [T] the signal is 
initially monochromatic during the growth phase of the 
mode and evolves successively to lower frequencies as the 
star spins down. We find that if an unstable /-mode 
with multipoles I = m = 3, 4 saturates at a sat = 10~ 4 , 
the gravitational-wave signal can be detected by ET for 
the most part of the instability evolution for a supra- 
massive star located in the Virgo Cluster at a distance 
of 20 Mpc. For the standard N = 1 model an unstable 
I = in = 4 /-mode may generate a detectable signal only 
if the star rotates near the Kepler limit O > 0.98 Ok- 

A higher magnetic field may affect significantly the 
characteristic strain of the /-mode, as the number of 
accumulated cycles N cyc decreases. The results for the 
Bp — 10 12 G case are shown in Figure [7] for both the 
N = 1 and N — 2/3 models. Although the gravitational- 
wave signal is slightly weaker than in the previous B v = 
10 11 G case, it is still detectable by ET. The effect of 
the magnetic breaking is particularly evident if the star 
gets unstable at lower rotation rates. For instance, the 
characteristic strain of the N = 1 model with an ini- 
tial rotation rate of O = 0.97 Ok is significantly smaller 
than for the low magnetised model. We have also cal- 
culated the characteristic strain for models with a mag- 
netic field strength of B p = 10 13 G and found that the 
gravitational-wave signal drops below the detectors' sen- 
sitivity curves. The results presented in this Section can 

be easily rescaled with the /-mode saturation amplitude 

7 1/2 
as h c ~ a s ' at . 

If in addition a r-mode is present as well, the results for 
the gravitational wave strain depend on the value of the 
saturation amplitude for this auxiliary mode; see Figure [8] 
for the N — 2/3 polytrope. Supposing a small r-mode 
saturation, the expected signal can still be detected with 
ET from a source in the Virgo cluster. However, as the 
saturation amplitude of the r-mode is increased, its insta- 
bility dominates the evolution of the /-mode instability. 
The star loses angular momentum too fast for the /-mode 
to grow substantially before it leaves its instability win- 
dow (see also Figure [6]), leading to a significant reduction 
in the gravitational wave strain for large r-mode satura- 
tion values. In these cases the signal of the /-mode will 
drop even below the ET sensitivity curve. However, esti- 
mating the gravitational- wave growth time of the r-mode 
with a uniform density stellar model |53j , we find that an 
unstable I = m — 2 r-mode can generate a gravitational- 
wave signal detectable by ET. This means that when 
both the /- and r-modes are simultaneously excited and 
the r-mode saturation amplitude is a s r at < 10~ 7 , ET 
should be able to detect the gravitational-wave signal 
emitted by both these two modes. Instead, for larger 
a s r at the r-mode will be still detectable but the /-mode 
signal should be below the ET sensitivity curve. 



; 1 1 1 1 1 1 1 1 
- N = 2/3 


| 1 | 1 | , | 1 . 


^ OID. K = 0.98 


Adv. LIGO 


r d = 20Mpc - 

_ 




- 


IT •- - 








a = 10" 5 

r 






a = 10 -7 

r 








— a-10- 6 " 

r 


' 








I I I I I I I I I I I I I I I I I I I 

0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 

v [ kHz ] 



FIG. 8: The characteristic strain generated by the I = m = 4 
/-mode instability when an / = m = 2 r-mode is excited 
as well. The star is a relativistic polytrope with N = 2/3, 
B p = 10 11 G and is located at 20 Mpc. It enters the r-mode 
instability window at Q — 0.98S7k- The Figure shows three 
cases with different r-mode saturation amplitudes while the 
saturation amplitude of the /-mode is fixed to a aat = 10 -4 . 



In Figure [7] (right panel) the characteristic strain fea- 
tures some vertical spikes. They arise from the emission 
of a quasi-monochromatic signal during the initial mode's 
growth phase and the use of equation ( |37[ ) to determine 
the accumulated mode's cycles in all the evolution. As 
already discussed in these spikes can be considered 
unphysical and we neglect their effects in our discussion. 

VI. CONCLUSIONS 

We have presented in this work the first dynamical 
study of the /-mode instability in General Relativity 
which incorporates the effect of viscosity, magnetic fields 
and unstable r-modes. These are the most dominant 
effects which may have a significant impact on the evo- 
lution of this mode. Our neutron star's models are rela- 
tivistic and described by a polytropic equation of state. 
We consider more in detail two sequences of uniformly ro- 
tating stars which can rotate up to the Kepler limit and 
describe, respectively, a neutron star with Mb = 1.4M Q 
and a supramassive neutron star with Mb = 1.6M . 

The /-mode instability may develop in the aftermath 
of a supernova explosion, when a new born neutron star 
may spin very rapidly. The most unstable /-mode is ex- 
pected in the I = m = 3 and 4 multipoles, which have 
a relatively large instability window. Considering var- 
ious cases in which a neutron star enters the instabil- 
ity window at different rotation rates, we find that the 
gravitational-wave signal emitted during the instability 
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may be detected by ET from a source located in the 
Virgo cluster, and as expected the gravitational signal 
is stronger for the more massive model with N = 2/3. 
This result is valid for an /-mode which saturates at 
E ~ 10~ 6 Mqc 2 and for a signal integration time < lyr. 
From the Virgo cluster we should expect about 30-60 su- 
pernova explosions per year, by assuming a rate of 2-3 
events per century in our Galaxy. This means that if 
few of them leaves behind a very rapidly spinning proto- 
neutron star, we might be able to detect these events 
from the gravitational radiation emitted during the /- 
mode instability. The number of these potential sources 
may be even more promising if we note that in more 
massive stars the gravitational-wave signal remains de- 
tectable for about 10 years. 

Another important result of our simulations is that 
the heat generated by shear viscosity during the /-mode 
saturation phase prevents the star from entering the 
regime of mutual friction. In fact, the shear viscosity 
re-heating balances the neutrino cooling and leads to a 
nearly isothermal evolution. The star therefore leaves 
the instability window at lower rotation rates, and even 
a moderate change in the superfluid transition tempera- 
ture does not change this result. 

The magnetic torque as well as an unstable r-mode 
may accelerate the transition of the star through the 
instability window of the /-mode. This may limit the 
growth of the mode's amplitude and therefore the strain 
of the gravitational-wave radiation. Our results show 
that the magnetic field affects the /-mode instability 
when B p > 10 12 G, and its influence is more relevant 
in the standard neutron star model. For the r-mode we 
find that it must reach the maximum value expected from 
non-linear mode coupling studies in order to affect con- 
siderably the /-mode evolution. However, a definitive 
answer to this issue cannot be given as the maximum 
/-mode amplitude is still unknown. 

In this study, we have not addressed the 'classical' I = 
ttl — 2 /-mode instability, which is potentially the lowest 
order unstable mode driven by gravitational wave emis- 
sion. As it was shown previously in [17], the prospects 
for a dynamical long-term evolution of a quadrupolar in- 
stability are rather pessimistic given the diminished size 
of the corresponding instability window. However, drop- 
ping the Cowling-approximation which we used here will 
certainly improve the situation for the I = m = 2 case 
since in full general relativity the critical angular veloc- 
ity for the onset of the CFS-instability is shifted towards 
lower values [4]. On the other hand, including relativis- 
tic effects for the r-mode might have the opposite effect 
in the sense that relativistic perturbation theory hints 
towards a weakening, i.e. much larger growth times, of 
the instability when compared to Newtonian results [54] . 
This also works in favor of an enhanced detectability even 
if /- and r-mode are unstable at the same time. 

So far we used polytropic equations of state for mod- 
elling relativistic neutron stars. A natural extension of 
this work should also incorporate a wide range of realistic 



equations of state. This most probably will not change 
the main results of this work but for sure will lead to 
alterations in mode frequencies, growth times and grav- 
itational wave detectability. At the same time, in more 
compact neutron star models one needs also to take into 
account the effect of direct Urea processes. These /3- 
reaction processes are energetically favourable when the 
proton fraction is sufficiently large, x p > 1/9, [55] . and 
they lead to a stronger bulk viscosity which may limit 
considerably the instability window of the /-mode [40] . 

Another potentially significant damping mechanism 
neglected here is the formation of a crust and its impor- 
tance on the /-mode instability. Additional dissipation 
will occur at a viscous boundary layer at the crust-core 
interface or more directly, the crust could break by a large 
amplitude /-mode. This would release a non-negligible 
amount of energy for dissipation at the fracture sites. 
Of course this mechanism depends on how fast a crust 
is formed after the creation of a hot, rapidly rotating 
proto-neutron star. Large amplitude oscillations might 
even delay the formation of the crust in the nascent neu- 
tron star. 
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Appendix A: The evolution equations of the /- and 
the r-mode instability 

In this Section, we derive the evolution equations for 
studying the evolution of the gravitational wave instabil- 
ity driven by two modes. Although these equations have 
been applied in Sec. |IV Cl for studying the /- and r- modes 
the formalism is general. As in the single mode case the 
main equations are ^ and (14 1, which now have to be 
extended to describe the instability of two modes. 

The first two equations can be derived from the mode 
energy variation ^ 



~dT 



din dn 

: dn ~dl 



2q x 



(Al) 



where the index x = f , r with x ^ y identifies the quan- 
tities related to the /- and r-mode, respectively, and 
we have defined the energy of each x-mode as E x = 
a x E x (Q). 

A third equation can be determined from the angular 
momentum evolution equation (14 1, with the canonical 
angular momentum now given by 



(A2) 
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where J* (f2) and J r c (fi) describe the dependence of the 
canonical angular momentum of each mode on the angu- 
lar velocity of the star. The other term that changes in 
equation ( 14 1 is the gravitational radiation torque, which 



now can be decomposed as 



dt 



2j[ 

V 

gw 



2JI 



gw 



(A3) 



where Tg W and t* w are the gravitational growth time of 
the /- and r-modes. Taking the time derivatives, the 
angular momentum conservation equation reads 



T t dctf 
Jc dt ' 

_2Jl 
r f 

• gw 



.da T 
~dt 



2j: dj, 



dJ s 

mag 



dn 



dJi 

' dQ, 



dn 

~dt 



1 gw 



dt 



(A4) 



As final step, we can combine equations ( Al I and ( A4) 



and obtain the following system of evolution equations 



da x 



2ox 

1 gw 

A, 



2a x D x 



2a x a Y A X F V 



D 



D 



dft 
~dt 



D 

E 2F x a x 

X '« 



dj s \ 1 dj,, 
~dn) ~~dt 



mag 



dJ s 
~dtt 



1 dJ, 



mag 



dt 



(A5) 
(A6) 



In equations ( A5 ) and ( A6 ) we have defined the following 
functions: 



.4 



d\nE* 



dVL ' 
Pyi — A x • F x ? 



F — T* 



dJ s 
~dtt 



Qx = 



dJ* fdJ 



dn \ dn 



(A7) 
(A8) 



where 



D = 1 + J2^(Q X -P X ) , (A9) 

X 

D X = D + a x P x . (A10) 



Similar to the single mode instability, we may deter- 
mine the evolution equations of a star in case one of the 
modes saturates. When the y-mode saturates its am- 
plitud e is nearl y co nstant, i.e. da y /dt = 0, and equa- 
tions ( A5 ) and ( A6 1 now read 



da x 

~dT 



D, 



dQ 
~dt 



2ox _ 

r x 

gw 
a x A x 

Dy 
I gw 



y 

-1 



D v 



2a x Dg + 2a x a y A x F y 

'v ^ 

dJ s 

~dn 

2F x a 



dJ n 



dt 



(All) 



+ 



dJ s 
~dtt 



dJ n 



dt 



(A12) 

where D q = 1 + J2 X «xQx- 

When both modes have reached their respective satu- 
ration amplitudes we approximately have that dat/dt = 
da T /dt — 0. As a result the star rotation evolves accord- 
ing to the following equation 



D, 



dn 

~dt 



EF x a x 
T x 

x 9W 



dJ s 
dVL 



dJ,, 



mag 



dt 



(A13) 



where in a low magnetised star both t and T r gw govern 
the spindown-rate of the star. 
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